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ABSTRACT 

Context. The formation of subdwarf B (sdB) stars is not well understood within the current framework of stellar single and binary 
evolution. 

Aims. In this study, we focus on the formation and evolution of the pulsating sdB star in the very short-period eclipsing binary 
PG 1336-018. We aim at refining the formation scenario of this unique system, so that it can be confronted with observations. 
Methods. We probe the stellar structure of the progenitors of sdB stars in short-period binaries using detailed stellar evolution cal- 
culations. Applying this to PG 1336-018 we reconstruct the common-envelope phase during which the sdB star was formed. The 
results are interpreted in terms of the standard common-envelope formalism (the a-formalism) based on the energy equation, and an 
alternative description (the y-formalism) using the angular momentum equation. 

Results. We find that if the common-envelope evolution is described by the a-formalism, the sdB progenitor most likely experienced 
a helium flash. We then expect the sdB mass to be between 0.39 and 0.48 M Q , and the sdB progenitor initial mass to be below ~2 
M Q . However, the results for the y-formalism are less restrictive, and a broader sdB mass range (0.3 - 0.8 M G ) is possible in this case. 
Future seismic mass determination will give strong constraints on the formation of PG 1336-018 and, in particular, on the CE phase. 

Keywords, subdwarfs - stars: evolution - stars: individual: PG 1336-018- binaries: close - binaries: eclipsing - methods: numerical 
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Subdwarf B (sdB) stars are the dominant population of faint 
blue objects at high galactic latitudes ( Green et al.| |T986), and 
are found in both the disk and halo. They are also ubiquitous 
in giant elliptical galaxies, where they are believed to be the 
main source of the ultraviolet excess ( |Brown et al.| 1997"] ). In the 
Hertzsprung-Russell diagram they lie on the blue extension of 
the Horizontal Branch, and are therefore also known as Extreme 
Horizontal Branch (EHB) stars. It is generally thought that they 
are low mass (0.5 M ) core-helium burning stars with extremely 
thin h ydrogen envelopes (< 0.02 Mg) (He ber|1986[|Saffer etaL 
|1994| l. Their envelopes are too thin to sustain hydrogen burn- 
ing, hence they will evolve directly to the white dwarf cooling 
track after core-helium exhaustion, without going through the 
Asymptotic Giant Branch and Planetary Nebulae phases. 

It is not clearly understood how the sdB progenitor manages 
to loose almost its entire hydrogen-envelope, but nevertheless 
starts core-helium fusion. Both single star evolution with en- 
hanced mass loss on the Red Giant Br anch (RGB) (|D'Cruz et aT] 
1996[ l, and binary evolution models (Men gel et al.||1976 l have 



been proposed as formation channels. Extensive surveys show 
that a large fraction of sdB stars are in binaries (e.g.|Allard et al 



|1994[ |Morales-Rueda et al.||20fjo) . This motivated" |Han et al. 



(2002 2003 i to perform a detailed investigation of the main 



binary evolution channels that can produce an sdB star. They 
found that an sdB star can be formed after one or two common- 
envelope (CE) phases producing a short-period binary (P = ~0. 1 
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to ~10 d) with respectively a main-sequence (MS) or a white 
dwarf (WD) companion. Only one phase of stable Roche Lobe 
overflow (RLOF) is predicted to contribute to the sdB popula- 
tion. This channel produces a wider binary (P — ~1 to ~500 d) 
with an MS companion. Single sdB stars are explained by the 
merger of two helium white dwarfs (WD). The binary popula- 
tion synthesis models for these formation channels ( |Han et aT] 
2003 ) predict a mass distribution of sdB stars that sharply peaks 
at the canonical value of 0.46 M , but it is much wider (0.3 - 
0.8 Mo) than previously assumed. The wide mass range is due 
to stars which ignite helium under non-degenerate conditions. 
These systems had not been explored as sdB progenitors before. 

A fraction of sdB stars show multimode short-period oscil- 
lations with amplitudes in the milli-magnitude range. They are 
observed to have surface gravities (logg) between 5.2 and 6.2 
and effective temperatures (T e ff) between 28,000 and 36,000 K. 
This class of pulsators is known as sdBV or V361 Hya stars. 
They are also often referred to as EC 14026 stars after the proto- 
type, discovered by Kil kenny et al.| ( [T997] >. Independently, these 
oscillations were theoretically predicted to be driven by an opac- 
ity mechanism (Charpi net et al.|1996[ ). A seismic study can pro- 
vide detailed constraints on the sdB interior, most importantly 
the total mass and the mass of the hydrogen-envelope, which are 
essential ingredients to tune the sdB formation scenarios. 

An excellent laboratory for a detailed seismic and evolution- 
ary study is the sdB pulsator in the short-period (2.4 h) eclipsing 
binary PG 1336-018, also known as NY Vir. The sdB primary 
was discovered to pulsate by |Kilkenny et al. ( 1 998|l, and has been 
the target of a Whole Earth Telescope campaign ( |Kilkenny et al.| 
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|2003| >. However, an adequate seismic model has not been deter- 
mined yet due to the lack of colour information. In future work, 
we will attempt to achieve this by using high-precision VLT pho- 
tometry and spectroscopy of this target star. An overview of the 
data and the orbit solution can be found in Vuckovic et al. (2 007) >. 
Here we present a progenitor study of PG 1336-018 from a the- 
oretical point of view. Assuming the sdB mass to be the canonical 
0.5 M , Kilke nny et al.| ([T998 ) derived a mass for the companion 
of 0.15 M Q and estimated its class to ~M5V. In view of the wide 



sdB mass distribution predicted by Han et al. ( 2003 i, we drop the 
assumption on the sdB mass and investigate the range of initial 
system parameters for binaries that evolve into a PG 1336—018- 
like configuration. 

The current orbital separation, ~0.8 R (Vuckovic et al. 
|2007| l, is much smaller than the radius the sdB progenitor had 
as a red giant. This implies that the system evolved through a 
common-envelope (CE) and spiral-in phase. At the start of mass 
transfer, the giant must have achieved a certain minimum core 
mass, for the core to still ignite helium after loss of the enve- 
lope. However, mass transfer must have started before the giant 
reached the tip of the RGB, because after the tip the giant starts 
to contract. The range that the giant's core mass and radius can 
have, has been calculated by Han et al. ( 2002| l as a function of 
the zero-age main-sequence (ZAMS) mass. They showed that 
the minimum core mass required for helium ignition is typically 
within 5% of the core mass at the tip of the RGB, where their 
definition of the core mass boundary is closely related to the 
layer of maximum energy production rate ( |Han et al.|1994 1. We 
perform a similar study here, but with some refined constraints. 
Most importantly, we take into account that the minimum core 
mass for helium ignition depends sensitively on the hydrogen 
envelope that is kept by the star. 

While a progenitor study of PG 1336-018 sheds light on 
the origin of sdB stars, the future evolution of this system 
is interesting in the context of cataclysmic variables (CVs). 
ISchreiber & Gansicke| ( [2003l calculated that PG 1336-018 will 
evolve into a semi-detached configuration within the Hubble- 
time and thus is representative for progenitors of present-day 
CVs. However, they mistakenly took the relative radius of the 
secondary from Kilkenny et al. ( 1998| l as the absolute radius. 
Furthermore, we now have more accurately determined system 
parameters than Kilkenny et al. ( 1998 1. Therefore we reinves- 
tigate the status of PG 1336-018 as a pre-CV as well and find 
that the main conclusion by Schreiber & Gansicke (2003) re- 
mains true, i.e. PG 1336-018 is representative for progenitors of 
present-day CVs. 

The outline of this paper is as follows. Section [2]briefly de- 
scribes the stellar evolution code used in this study, and the pro- 
cedure we adopted for our calculations. Section [3] presents the 
results we obtained by probing the stellar structure of progen- 
itors of sdB stars in short-period binaries at the onset of mass 
transfer. In particular, we present the pre-CE orbital separation 
of possible progenitors of PG 1336-018. This is used in Section 
13.21 to constrain the CE evolution. The results are discussed in 
Section|4]and summarized in Section|5] 



2. The evolutionary calculations 

2. 1 . The stellar evolution code 

We compute the stellar evolution with the numerical computer 
code originally developed by Eggleton] ([197 1 



Eggleton 'etaLl ( |1973| l and updated by |Han et al.| ( |1994| l andfPo" 



1972 1973 



et al.| ( |1995| [1998|J The updated version of the code uses an 
equation of state that includes pressure ionization and Coulomb 
interaction, opacity tables derived from |Rogers & Iglesias ( 1992 ) 

994), nuclear reactio n rates from 



and |Alexander & Ferguson| 



Caughlan et al. (1985 1 and 



Caughlan & Fowler (1988), and 



neutrino loss rates from |Itoh et al.| ( |1989| |1992) . The code 
uses a self-adaptive, non-Lagrangian mesh. During an itera- 
tion it solves implicitly and simultaneously the stellar structure 
equations, the chemical composition equations and the equa- 
tions governing the mesh-spacing. Both convective and semi- 
convective mixing are treated as diffusion processes. Izzard & 



Glebbeek (2006) recently developed a graphical user interface, 
Window To The Stars (WTTS), to Eggleton's code. WTTS sig- 
nificantly simplifies running the code and allows immediate 
analysis of results. 

We use a mixing-length parameter (the ratio of the mixing- 
length to the local pressure scaleheight) of a = l/H p = 2.0. 
Convective overshooting is included using an overshooting pa- 
rameter 6 ow =0.12 which corresponds to an overshooting length 



of ~ Q.25H,,. We use a Reimers' wind mass-loss rate (|Reimers 



[1975]), 



= 4x 10 



. 13 (R/R Q )(L/ L q ) _ 
71 (M/M B ) [M ° yr 



1, 



(1) 



with an efficiency of rj = 0.4 ( |Iben & Renzini 1983 ; Carra roet al.| 
|1996[ ). The metal l icity i s taken to be Z = 0.02. 

In Han et al. ( 1994[ ), the stellar core boundary is related to 
the layer of maximum energy production rate. This will not be 
applicable for EHB stars since they have very thin inert hydrogen 
envelopes. Therefore, we define the core to be the inner region 
with a hydrogen mass fraction X < 0.10. This will generally 
lead to lower values for the core mass, although it will not give a 



significant difference for degenerate cores (Dewi & Tauris 2000 



Tauris & Dewi 2001 ). In principle, it is not important how the 
core boundary is defined, provided that it not assumed to be the 
bifurcation point above which all the material is ejected. 

2.2. Procedure 

We have used the Eggleton code to follow the evolution along 
the RGB of stars with ZAMS masses in the range 1-4 M . 
We do not study more massive progenitors, because they will 
result in EHB stars of mass above 0.8 M Q Pan et al.[|2002| ), 
which are too hot to become sdB pulsators. We approximate 
the CE phase by removing the envelope at a rate of 10~ 6 M* 
yr~' where M* is the mass of the star, while keeping the com- 
position constant. Clearly, this is a crude approximation, as CE 
evolution involves much higher rates of mass loss causing the 
star to lose hydrostatic equilibrium. It is, however, expected that 
the subsequent evolution does not depend on the mass loss his- 
tory, but it is mainly determined by the amount of hydrogen left. 
Unfortunately, it is unknown how much envelope is ejected dur- 
ing the CE phase. What we do here is derive an upper limit for 
this quantity. First of all, we expect the star to contract after CE 
ejection. However, as sdB spectra are generally dominated by 
hydrogen lines, we do not remove more envelope than down to 
X = 0.10. Secondly, as we are investigating post-CE sdB stars, 
the stars cannot have a too large radius after mass-loss, therefore 
we require /? pos t-CE < 10 2 R . This is still a much too large limit 
for PG 1336-018, but it might apply to other post-CE sdB sys- 
tems, and for now we want to keep the discussion general. The 



1 A write-up of the most recent version of this code can be obtained 
from P. Eggleton at ppe@igpp.ucllnl.org. 
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final condition we impose is that the star must evolve to the EHB 
with temperatures above 28,000 K, as is characteristic for pul- 
sating sdB stars. Which of these three conditions will determine 
the maximum remaining envelope depends on the situation. To 
clarify this, we state the four different scenarios with the domi- 
nating criterion: 

1) M Z ams < 2 M Q , near M core 

,min • 

When the core is degenerate, the core boundary is very dis- 
tinct. In this case, it is reasonable to suppose that nearly all 
the material above the very compact core is expelled. Also, 
the large surface gravity at the core boundary will restrict 
the remaining hydrogen envelope to be M env < ~10~ 3 M Q , 
for exact values see Table[T] We found that thicker envelopes 
will continue to burn hydrogen and these models have too 
large radii (~150 Rq) to fit in the narrow orbits of post-CE 
systems. 

2) M Z ams < 2 Mq, near RGB tip: 

When the core degeneracy is liftecj^] during the CE ejection, 
we find that the remaining envelope can have a mass up to 
~10~ 2 M Q (see Table [T) , consistent with the evolutionary 
studies of |Caloi| ( fT989| T and |Dorman et al.| ( fT993] ). In this 
case, we determined the maximum envelope which allows 
the star to reach the EHB with an effective temperature above 
28,000 K. 

3) M ZAMS > 2 M , near RGB tip: 

When these stars with non-degenerate cores are close to he- 
lium ignition, the situation is comparable with case 2. Also 
here, the maximum envelope mass (~10- 2 M , see Table[l) 
follows from the requirement that the sdB star must reach 
T e g > 28, 000 K during core helium burning. 

4) M Z ams > 2.5 Mq, near M core , min : 

These stars can ignite helium already when they loose their 
envelopes at the end of the MS. For the end-of-MS models 
we have not been able to derive a realistic upper limit for the 
remaining hydrogen envelope, because they tend to expand 
even when we remove the entire envelope, i.e. down to X — 
0.10. For these models, we assumed that the entire envelope 
was ejected. 

Using the above conditions for how much envelope should 
at least be removed, we determined the minimum core mass for 
helium ignition, and the stellar structure at onset of mass transfer. 

3. Results 

3.1. At the onset of mass transfer 
3.1 .1 . The stellar structure 

We are interested in the stellar structure of the sdB progenitors at 
the onset of mass transfer. In particular the total mass, the radius, 
the remnant mass, and the binding energy of the removed enve- 
lope are important for constraining the CE evolution. In Table [T] 
we present the relevant stellar parameters corresponding to the 
minimum core mass for helium ignition and the tip of the RGB. 
We also give limits on the total mass and the mass of the hy- 
drogen envelope of the post-CE star, which will be close to the 
values of the sdB star itself. 

Note that for Mzams < 2 M Q , the remaining hydrogen en- 
velope can be thicker when mass transfer started at the tip of 
the RGB. In this case we removed the envelope of post-He- 
flash models, i.e. we assumed that the core degeneracy was lifted 
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Fig. 1. Evolutionary tracks of sdB progenitors until the tip of 
the RGB in the HR diagram. At the upper two tracks the red 
giant ignites helium quiescently, while at the lower two tracks 
the helium flash occurs. On the dotted part of the track, the core 
mass is above the minimum required for helium ignition. Thus 
when the giant looses its envelope in this stage of the evolution, 
it may become an sdB star. 



during the CE ejection. For the models at the end of the main- 
sequence, i.e. Mzams > 2.5 M Q and M core = M core m i n , we could 
not determine the hydrogen envelope reliably (see Section pT4) . 
Instead we removed the envelope down to X = 0.10. 

It is interesting to note that the stellar structure at the onset 
of mass transfer is quite different depending on whether the gi- 
ant experienced a helium flash or ignited helium quiescently. For 
degenerate cores, core contraction is inhibited by the degeneracy 
pressure, thus mass transfer can only have started very near the 
tip of the RGB, for the core to still ignite helium. More mas- 
sive stars (Mzams > 2 M ) achieve temperatures high enough to 
avoid core degeneracy. They do not have to be close to the tip 
of the RGB at the onset of mass transfer to still ignite helium, 
as long as their core mass exceeds the absolute minimum for he- 
lium ignition (-0.3 MgQ. IHan et al.| ( |2002| l found that stars with 
Mzams ^ 2.5 M Q will burn helium even when the envelopes 
are lost when passing through the Hertzsprung gapj^This can be 
clearly seen in the HR diagram shown in Fig.[T] 



3.1.2. The orbital separation 

It is generally assumed that mass transfer starts as soon as the 
giant fills its Roche lobe, i.e. when its radius equals the Roche 
radius. The Roche radius is approximated by 



Rl = 



QA9q 



2/3. 



0.6q 2 ' 3 +ln(l + q 1 ' 3 )' 



(2) 



where q — M1/M2 is the mass ratio and a the orbital separa- 
tion (Eggleton 198 3). Howe ver, if the binary system is tidally 
unstable ( Counselman| 1 973) , the Roche geometry is not appli- 
cable. It is expected that the time-scale of tidal evolution is much 
shorter than the time-scale of nuclear expansion, i.e. when a tidal 
instability sets in, the stars spiral inwards before the stellar struc- 
ture can change significantly. Therefore, when a tidal instability 



2 We used approximated post-He-flash models, as the code cannot 
calculate through the helium flash (Pols et al. 1998 1. 



Stars with 2 < Mzams < 2.5 M Q leave the main-sequence with core 
masses 

they are quite near the RGB tip. 
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Table 1. The stellar structure before and after mass transferF] 
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" For each ZAMS mass, the first line gives the stellar structure corresponding to the minimum core mass for helium ignition. The second line 
corresponds to the tip of the RGB. The columns are respectively: Mzams = zero-age main-sequence mass; M, = total mass of giant; M core = helium 
core mass of giant, R, = radius of giant; A. p and A tat are dimensionless parameters indicating respectively the gravitational binding energy and the 
total (including thermal) binding energy of the ejected envelope (see Section |3.2.1fr ; / = the moment of inertia; M cnv = hydrogen envelope left after 
CE ejection; M remnlmt = remnant mass after CE ejection. 



sets in before the giant fills its Roche lobe, this can also cause a 
CE. The binary becomes tidally unstable when the spin angular 
momentum of the stars exceeds one third of the orbital angular 
momentum h: 



(h + h)tj > -h 



with 



GaM 2 Ml 

giant 2 



giant 



+ M 2 



(3) 



where I\ and I 2 are the moments of inertia of the stars and a> is 
the angular velocity ( Hut|1980 1. Since the giant is much larger 
and more massive than the MS companion, we have /; » I 2 . 
The moment of inertia of the giant, I\, follows from detailed 
evolutionary calculations, and can be found in Table[T] 

Using the allowed range for the radius and mass of the pri- 
mary at the onset of mass transfer (Table [T}, we can calculate 
the pre-CE orbital separation a, for a given secondary mass. If 
the binary is tidally stable when it fills its Roche lobe, Eq. (|2]i 
can be used to obtain a, . Otherwise, the CE is formed by a tidal 
instability and we use 



I\U) 



1 



(4) 



In Fig. [2] the range of a, is given as a function of the primary 
ZAMS mass. We have chosen a secondary mass of 0.12 M Q 
here, which is a best-fit value for PG 1336-018, observation- 
ally derived by Vuckovic et al. ( |2007[ ). The effect of a different 
secondary mass is investigated in section |4~T| 



350 
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"q 200 
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a 150 
100 
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He flash 



Fig. 2. The pre-CE orbital separation a, as a function of the pri- 
mary ZAMS mass for M 2 = 0.12 M G . At the dotted line, mass 
transfer started when the helium core mass reached the mini- 
mum required for helium ignition. At the solid line, mass trans- 
fer started at the tip of the RGB. Thus the shaded area gives the 
possible values of a, for PG 1336-018. At the triangles, the CE 
is formed by dynamically unstable RLOF. At the squares, the CE 
is caused by a spiral-in due to a tidal instability. 



3.2. Common-envelope ejection 

3.2.1. The energy equation: the a-formalism 

In the original spiral-in picture ( |Paczynskr]|1976 1, the compan- 
ion experiences drag forces as it moves into the envelope of the 



giant and spirals inwards. This idea is applicable when a tidal 
instability causes the CE or sets in soon after the CE is formed, 
i.e. when the mass ratio is high. In this case, the CE is not in 
co-rotation and the envelope is heated and unbound by friction. 
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The orbital energy E OI b released in the spiral-in process, is used 
to eject the envelope with an efficiency a, 



^(^arb,f ^orb,i) — ^ei 



(5) 



where £ env is the binding energy of the ejected envelope, and 
the subscripts i and / denote values before and after the CE 
phase. In principle, one expects < a < 1 . However, in order to 
explain observed binaries one often finds that a exceeds unity. 
This might indicate that other energy sources can contribute to 
the ejection of the envelope, e.g. the luminosity of the giant. 
Still, a very high value for a is not anticipated, since it would 
be physically difficult to explain where such a large amount of 
energy could come from. The poorly understood physics of the 
CE phase does not allow us to set a hard limit on a, but here we 
assume < a < ~5 to be realistic. 

It is reasonable to suppose that the secondary did not accrete 
any matter since the mass transfer time-scale is short. Expression 
|5]) can then be written as 



-GM remnant M 2 GM giant M 2 



2af 



?giantM 2 \ _ 

2a\ I 



GM glant M elw 



AR 



(6) 



giant 



where we have expressed E env in terms of the structural parame- 
ter A (Webbink |1984| l. It is straightforward to calculate the com- 
bined parameter aA from Eq. (|6|. To isolate a one usually takes 
A - 0.5 ( |de Kool|T99 ), but an exact calculation should take into 
account that A depends on the stellar structure. 

The total binding energy consists of the gravitational binding 
energy and the internal thermodynamic energy U, 



r M ^' i gm \ 

+ U)dm. 



(7) 



Taking the total binding energy of the envelope to calculate A 
implies that the entire internal energy is used efficiently in the 
ejection process. However, it is uncertain how much of the in- 
ternal energy contributes to the ejection of the envelope. This 
uncertainty is expressed in a parameter a±, introduced by |Han 
|etaL| ( fT99^ , 



■JM, cmmml V 



GM\ 



dm + a t h 



Jm,„ 



Udm. 



(8) 



Expression ([8]) can be regarded as the effective binding energy 
of the envelope and is used to derive A. We calculated A for 
a th = 1 and ar tn = 0, representing respectively the total bind- 
ing energy (A tot ) and the gravitational binding energy (A gr ), see 
Table [T] In Fig. [3] a is plotted as a function of the sdB mass for 
A = A„. In principle, we do not expect the ionization energy to 



contribute to the ejection of the CE, as |Harpaz] ( |1998 1 argued 
that after recombination the opacity drops sharply, hence the re- 
leased energy flows outward without unbinding the material. In 
the case that the internal energy of the envelope does contribute 
to the ejection process (i.e. or t h = and A = A gr ), this will result 
in a lower (~factor 2) value of a, as is expected from the virial 
theorem. This can be seen in Fig. |4] 

It is shown that, in general, non-degenerate helium ignition 
requires unphysically high a-values (a > 5). However, if the 
internal energy of the envelope is efficiently used to unbind 
the envelope, we find a > 2. Us ing the initial mass function 

1993} and the d istribution in orbital 
' " 1978) , we found that the 



<£(M) °c M 2 7 ( |Kroupa et al. 
separation T(a) oc a" 1 (|Kraicheva et al 



number of systems that experienced a helium flash is comparable 
to the number of systems that ignited helium non-degenerately 
with a < 5. Thus, if the internal energy can unbind the envelope, 
we are less confident in excluding the non-degenerate scenario. 
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Fig. 3. The CE parameter a as a function of the sdB mass for 
M2 - 0.12 Mq. At the dotted line, mass transfer started when the 
helium core reached the minimum required for helium ignition. 
At the solid line, mass transfer started at the tip of the RGB. 
The shaded region in between indicated the possible a-values for 
PG 1336-018. We assumed here that the binding energy of the 
envelope is determined by the gravitational energy, i.e. A = A^. 
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Fig. 4. The CE parameter a as a function of the sdB mass for 
M2 = 0.12 M . At the upper solid line, we have A = /l gr and at 
the lower dashed line, we have A = A lot . We assumed that mass 
transfer started when the red giant reached the tip of the RGB. 



3.2.2. The angular momentum equation: the y-formalism 

Nelemans et al. (2000, 2005) found that the first phase of mass 
transfer of observed double white dwarfs cannot be described 
by the standard a formalism, nor by stable RLOF. They argued 
that, for binaries with mass ratio close to unity, the common en- 
velope is formed by a runaway mass transfer rather than a decay 
of the orbit. In this case, the angular momentum of the orbit is 
so large that the common envelope is brought into co-rotation. 
Consequently, there are no drag forces that can convert orbital 
energy into heat and kinetic energy. |Nelemans et al. ( 2000| l de- 
scribed this scenario in terms of the angular momentum balance, 
the y-formalism. The assertion is that the specific orbital angu- 
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lar momentum carried away by the envelope is y times the initial 
specific orbital angular momentum, 



Ji 



M, 



giant 



+ M 2 



Ji-Jf 



(9) 



Although the -/-formalism was originally developed for double 
white dwarfs, Nelemans & Tout (2005) showed that most ob- 
served sdB binaries can also be explained by y ~ 1 .5. The phys- 
ical motivation for this empirical description might be super- 
Eddington mass transfer ( Beer et al.|2007| l. Originally this idea 
was put forward to explain systems in which a main-sequence 
star transfers mass to a neutron star or black hole ( jKing &| 
Begelman 1999}. Beer et al. ( 2001) proposed that this physical 
picture is also applicable to systems where a red giant overflows 
onto a main-sequence star. Their treatment gives an upper limit 
for y, since the specific angular momentum carried away by the 
envelope cannot be higher than the specific angular momentum 
of the secondary, 



ymax 



_ M giant + M 2 



M, 



env 



(M giant + M 2 ) 2 /_^env\ 

We have calculated the possible y- values for PG 1336-018 
and plotted them in Fig. [5] as a function of the sdB mass. We 
also plotted y max corresponding to the case that the envelope is 
ejected with the specific angular momentum of the secondary. 

The common value y ~ 1 .5 is found in the helium flash re- 
gion, but when the sdB mass approaches the canonical value of 
0.47 M G , the y-value increases steeply to > y max . This is be- 
cause the corresponding progenitor mass (Mzams = 1 M©) and 
the expelled envelope mass are low, thus the angular momen- 
tum carried away per unit mass is high. If the progenitor mass 
of PG 1336-018 is higher than 1 M Q , the y- values agree with 
the idea that the envelope is ejected with almost the specific an- 
gular momentum of the secondary. For non-degenerate helium 
ignition we find 1.1 < y < 1.2, and we cannot rule out this pos- 
sibility as an explanation for PG 1336-018. 



4. Discussion 

4. 1 . The influence of the secondary mass 

Assuming Newtonian mechanics, the radius of the secondary is 
given by 



4.2. Similar systems 

Up to now 39 sdBV stars are known, of which 16 are known 
binaries. 13 have spectroscopic F-G companions, two have in- 
visible WD companions, detectable only by radial velocity vari- 
ations of the sdB star, and only PG 1336-018 has been found 
to have an M-dwarf companion. Additionally, one star has been 
found to have a long period Jupiter mass planet (HS 2201+2610 
Silvotti et al.| |2007) detected only by the variation of the main 



pulsation period. Of the remaining 22 stars only four have been 
carefully checked for RV variations, the rest remains unexplored. 

Since PG 1336-018 is the primary of an eclipsing system, it 
is an exceptionally promising candidate for a study of the sdB 
binary formation channel. A few more non-pulsating sdB bi- 
nary systems with an M dwarf companion have been detected 
(see Table [2}, either through eclipses or strong reflection ef- 
fects. The eclipsing binaries HW Vir (Menz ies & Marang[l 986) 
and HS 0705+6700 ( |Drechsel ef"aT1|2001| ) show striking simi- 
larities with PG 1336-018. The fourth known eclipsing system, 
HS 223 1 +244 1 , while similar with respect to the short orbital pe- 
riod and other light-curve properties, turns out to have a very low 
mass M dwa rf companion, just at the substellar limit (Ostensen 
et al. 2007a)[j] A fifth very similar sdB+dM system was also re 



cently detected (Polubek et al.|2007|> in the OGLE-II photometry 
of the galactic bulge ( |Wozniak et al.|2002[ ), but no spectroscopic 
data has yet been obtained that can constrain the physical param- 
eters of this system. 

The very narrow distribution in orbital periods and compo- 
nent masses suggests that there is a common mechanism for 
producing such systems and that the CE phase of these sys- 
tems must have been very similar. To illustrate this, we have 
calculated a and y for the systems HW Vir, HS 0705+6700 and 
HS 2231+2441. For simplicity, we assumed here that the entire 
envelope was lost when the giant was at the tip of the RGB, and 
that the internal energy of the envelope cannot be used to unbind 
the envelope (i.e. a± = and A = A gI ). 

We can see in Table[2]that the CE parameters are indeed quite 
similar for these systems. The exception is HS 2231+2441, due 
to its very low mass companion. Note again that the canonical 
sdB mass (0.47 - 0.48 M Q ) corresponds to high y- values. 

Fig. [7] shows the position of PG 1336-018 in a - logg 
diagram, together with the other three systems of HW Vir type 
for which a spectroscopic temperature and gravity have been de- 
rived, among other sdB stars. Evolutionary tracks from Kawaler 
|& Hosfler|pb05) are also shown. We see that if HS 2231+2441 
as an sdB system, it is quite evolved on the EHB, while the other 
three systems are still in their first half of the EHB evolution. 



R 2 



R\r 2 



GM sdB r 2 
8 n' 



(11) 



The best-fit orbital solutions of Vuckovic et al. ( 2007 1 give av- 
eraged relative radii r\ = 0.19, r 2 = 0.21, and logg = 5.77. 
Substituting these values and the sdB mass range M sc jb ~ 0.3 - 
0.8 M Q in Eq. {□} gives R 2 ~ 0.13-0.21 R Q . Theoretical models 



of low-mass main-sequence stars by Baraffe et al. ( 1998 1 predict 
the mass range M 2 ~ 0.10 - 0.20 M Q for these radii. In Fig [6] 
the influence of the secondary mass on the a- and y-parameter 
is depicted. We see that a is lower when the secondary is more 
massive. However, the best-fit orbit solutions of Vuc kovic et al.| 
(2007 ) all have low masses for the secondary (M 2 — 0.1 1 - 0.13 
M ). The secondary mass has negligible influence on the y- 
parameter. 



4.3. The future evolution of PG 1336-018 



Following Schreiber & Gansicke (2003) but using updated sys- 
tem parameters, we determine the expected future evolution of 
PG 1336-018. The orbital period will decrease as the system 
loses angular momentum, 



./ 
1 



P 

3P' 



(12) 



For M 2 < 0.3 M , the secondary is fully convective and angular 



momentum loss via magnetic braking is negligible (Verbunt & 



|Zwaan|1981|). Thus we only need to consider angular momentum 



4 We note that this system might also be a non-He-burni ng post-RGB 
star with a substellar companion (0stensen et al. 2007b), but for our 
discussion we will assume it is an sdB system. 
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Fig. 5. The parameter y as a function of the sdB mass for Mi = 0.12 M . For visibility we have plotted the y- values for stars that 
ignited helium (non-)degenerately in the (right) left panel. At the lower dotted line, mass transfer started when the helium core mass 
reached the minimum required for helium ignition. At the lower solid line, mass transfer started at the tip of the RGB. The shaded 
region in between indicated the possible y-values for PG 1336-018. The upper dotted line corresponds to y max at the minimum core 
mass, and the upper solid line to y max at the RGB tip. 
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Fig. 6. The parameters a and y as a function of the sdB mass, plotted respectively in the left and right panel. The colour gradient 
gives the secondary mass. We assumed that mass transfer started at the tip of the RGB, and A - A gI . 



Table 2. System parameters of sdB+dM binaries similar to PG 1336-018.^] 



Ref. 



P(d) M sdB (M ) M 2 (M ) I hoj/h a 



HWVir 

PG 1336-018 



Wood & Safferl 19991 

VuckovTFeTand2U07b 
II T 



0.117 0.48 
0.101 0.39 
0.47 

0.53 

HS 0705+6700 |Drechsel et aL|(2001| 0.096 0.48 
HS2231+2441 |0stensen efaf] ( |200/a| > 0.111 0.47 



0.14 
0.11 
0.12 
0.13 
0.13 
0.075 



0.10 
0.80 
0.11 
0.17 
0.10 
0.20 



y 



0.34 2.81 

3.09 1.29 

0.36 2.76 

14.4 1.18 

0.32 2.79 

0.59 2.61 



" The a- and y-values given here correspond to the case that the entire hydrogen envelope was ejected at the tip of the RGB. Whenever the 
sdB mass allows it, we assumed that helium was ignited degenerately in a flash, indicated by normal font style. The CE parameters in boldface 
correspond to non-degenerate helium ignition. If l\ii)jh < 1/3, the CE is assumed to be formed by dynamically unstable RLOF. If l\U>lh > 1/3, 
the CE is assumed to be caused by a tidal instability. 
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Fig. 7. The location of PG 1336-018 in the r eff -log g plane 
plotted together with the other known eclipsing systems, non- 
eclipsing reflection binaries and some of the pulsating and non- 
pulsating sdB stars. The squares show known pulsators and the 
+ symbols non-pulsators. Also shown are evolutionary tracks for 
a sample of EHB stars (Kawaler & Hos tler|2005 1. 



loss due to gravitational radiation that follows from Einstein's 
quadmpole formula, 



32G 7/3 



2tt\ 7 / 3 



M\M\(M X +M 2 r 2/3 (y) 



(13) 



Eventually the secondary will fill its Roche lobe, initiating a sec- 
ond RLOF phase. As the secondary is essentially unevolved, we 
assume that its expansion due to its nuclear evolution is negligi- 
ble. Using Equation (|2l), Kepler's thir d law and the system pa- 

^ ( |2007| p] Mi = 0.466 M B , 
= 0.16 R Q , the orbital pe- 



rameters obtained by Vuckovic et al 
M 2 = 0.122 Mo, R\ = 0.15 R G , R 2 
riod of the semi-detached system can be calculated: P s d = 0.069 
d. The time from now until the secondary fills its Roche lobe 
follows from equations ( 12 1 and ( 13 1, 



fsd - 



5c 5 



(Mi + M 2 ) 



1/3 



256G 5 / 3 (2tt) 8 / 3 MiM 2 



-op; 



8/3 
now 



(14) 



We find ? s d = 9.9 x 10 s yr. Since the time-scale on which the sdB 
star evolves into a white dwarf is a few 10 8 yr, PG 1336-018 
will evolve into a short-period cataclysmic variable (CV). The 
present age of PG 1336-018 is maximally 12 Gyr. This is for 
Mzams = 1 M and the age decreases steeply with increasing 
mass. Thus the maximum total time until the system becomes 
semi-detached is less than the Hubble time. We conclude that 
PG 1336-018 is representative for progenitors of present-day 
CVs. Not surprisingly, we find that the similar systems given 
in Table|2]are also pre-CV candidates. 

4.4. Determining the minimum core mass and the remaining 
hydrogen envelope 

We cannot model the detailed physics of CE evolution as this is 
a hydrodynamical event. As an approximation, we used a low 
enough mass loss rate to ensure numerical stability, and pre- 
vent composition changes due to nuclear burning. However, in 



5 They found three statistically equivalent solutions, here we use their 
Model II. The results will not differ significantly for the other two mod- 
els. 



our models the composition profile does change slightly during 
the removal of the envelope. We are not sure if this is merely 
a numerical effect, or that some convective mixing might occur 
during CE ejection. In any case, we checked that resetting the 
chemical compositions to the values of the model before the CE 
phase does not change our main results for the a- and y-values. 
We plan to look into this issue in detail later, because it may be 
important for the seismic behavior of these stars. 

We have determined the rninimum core mass for helium ig- 
nition independent from Han et al. (2002 ), because we use a dif- 
ferent definition of the helium core mass. Moreover, we also are 
interested in how much of the hydrogen envelope is allowed to 
remain after CE ejection, so that we can make a consistent esti- 
mate of the binding energy of the ejected envelope. Our results 
are in general consistent with theirs, except for the minimum 
core mass for the helium flash. The reason is that we strip off 
more hydrogen, to exclude models that still have a growing core 
after CE ejection, because the star is then too large to fit in a nar- 
row post-CE orbit. Therefore, we find slightly higher minimum 
core masses for the helium flash than |Han et al.| (|2"002). 

Further, we observed that the models that can ignite helium 
when the envelope is lost during the Herztsprung gap, i.e. those 
with Mzams > 2.5 M , tend to expand even when we strip off the 
hydrogen envelope down to X = 0.10, because, at this stage, the 
hydrogen profile is rather smooth due to the shrinking convective 
core during the MS. So there is still some hydrogen left, mixed 
in the core, which continues to burn after CE ejection, causing 
the star to expand further. An interesting consequence is that the 
stars ejecting the CE at this earlier stage, might have more he- 
lium in their remaining envelope. We also plan to examine this 
possibility in detail in future work. 

5. Conclusions 

We have studied the stellar structure of the progenitors of sdB 
stars in short-period binaries. The narrow range a giant's core 
mass can have to still ignite helium after loss of its envelope, al- 
lows us to estimate the orbital separation at the onset of mass 
transfer. This enables us to constrain the CE phase. We have 
compared two different CE-descriptions; the a-formalism based 
on the energy equation (implicitly assuming angular momentum 
conservation), and alternatively, the y-formalism based on the 
angular momentum equation (implicitly assuming energy con- 
servation). Although we focused our study on the interesting sys- 
tem PG 1336-018, the methods we present here can readily be 
applied to other sdB stars formed in the CE ejection channel. In 
particular, Table[T]is very useful for such studies. 

Adopting the a-formalism implies that PG 1336-018 ignited 
helium in a degenerate flash. The sdB mass must then be be- 
tween 0.39 and 0.48 M Q . However, the results are less con- 
vincing if the internal energy contributed to the ejection of the 
envelope. Furthermore, the y-formalism does not rule out the 
possibility that the sdB progenitor ignited helium under non- 
degenerate conditions. In these cases, the possible mass range 
of the sdB star is wider: 0.3 - 0.8 Mq. If the CE is caused by a 
tidal instability, the a-formalism might be preferred. If the CE is 
formed by runaway mass transfer, the y-formalism is perhaps 
more likely ( |Nelemans~e t al. 2000). Unfortunately, a detailed 
physical description for the CE evolution is still missing and 
these conditions are not definite. For the time being, we consider 
both formalisms to be valid. 

We plan to use the pulsational properties of PG 1336-018 to 
make an independent high precision mass determination of the 
sdB star and its envelope. Clearly, this will shed more light on the 
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poorly understood CE phase. Also, if we are able to determine 
whether the sdB experienced a helium flash or ignited helium 
quiescently, this would provide us with much insight into the 
CE evolution. It is not clear however, if the helium flash leaves a 
(seismic) detectable imprint on the sdB interior. We will investi- 
gate these questions in a follow-up paper. 
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